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Popular Summary 

The global hydrological cycle is central to the Earth's climate system, with rainfall and the 
physics of precipitation formation acting as the key links in the cycle. Two-thirds of global 
rainfall occurs in the tropics with the associated latent heating (LH) accounting for three- 
fourths of the total heat energy available to the Earth's atmosphere. In the last decade, it has 
been established that standard products of LH from satellite measurements, particularly 
TRMM measurements, would be a valuable resource for scientific research and applications. 
Such products would enable new insights and investigations concerning the complexities of 
convection system life cycles, the diabatic heating controls and feedbacks related to meso- 
synoptic circulations and their forecasting, the relationship of tropical patterns of LH to the 
global circulation and climate, and strategies for improving cloud parameterizations in 
environmental prediction models. 

However, the LH and water vapor profile or budget (called the apparent moisture sink, or 
Q2) is closely related. This paper presented the development of an algorithm for retrieving 
Q2 using TRMM precipitation radar. Since there is no direct measurement of LH and Q2, 
the validation of algorithm usually applies a method called consistency check. Consistency 
checking involving Cloud Resolving Model (CRM)-generated LH and Q2 profiles and 
algorithm-reconstructed is a useful step in evaluating the performance of a given algorithm. 
In this process, the CRM simulation of a time-dependent precipitation process (multiple-day 
time series) is used to obtain the required input parameters for a given algorithm. The 
algorithm is then used to “ reconstruct ” the heating and moisture profiles that the CRM 
simulation originally produced, and finally both sets of conformal estimates (model and 
algorithm) are compared each other. 

The results indicate that discrepancies between the reconstructed and CRM-simulated 
profiles for Q2, especially at low levels, are larger than those for latent heat. Larger 
discrepancies in Q2 at low levels are due to moistening for non-precipitating region that 
algorithm cannot reconstruct. Nevertheless, the algorithm-reconstructed total Q2 profiles are 
in good agreement with the CRM-simulated ones. 
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Abstract 


The Spectral Latent Heating (SLH) algorithm was developed for the Tropical 
Rainfall Measuring Mission (TRMM) precipitation radar (PR) in part I of this study 
and improved in part II to estimate apparent heat source (Q x ) profiles. In this 
paper, the SLH algorithm is used to estimate apparent moisture sink (Q 2 ) pro- 
files. Procedure of Q 2 retrieval is the same as that of heating retrieval (Part I, II) 
except for using the Q 2 profile lookup tables derived from numerical simulations 
of tropical cloud systems from the Tropical Ocean Global Atmosphere (TOGA) 
Coupled Ocean-Atmosphere Response Experiment (COARE) utilizing a cloud- 
resolving model (CRM). 

The 0 2 profiles were reconstructed from CRM-simulated parameters [precip- 
itation top height (PTH), precipitation rates at the surface and melting level, and 
rain type] with the COARE table and compared then to CRM-simulated “true” 
Q 2 profiles, which were computed directly from the model water vapor equation. 
COARE, GATE, SCSMEX, and KWAJEX periods were used for the consistency 
check. The consistency check indicates that discrepancies between the SLH- 
reconstructed and CRM-simulated profiles for Q 2 , especially at low levels, are 
larger than those for Q x . Larger discrepancies in Q 2 at low levels are due to moist- 
ening for non-precipitating region that SLH cannot reconstruct. Nevertheless, the 
SLH-reconstructed total Q 2 profiles are in good agreement with the CRM-simulated 


ones. 



The SLH algorithm was applied to PR data, and the results compared to Q> 
profiles derived diagnostically from SCSMEX sounding data. Although discrepan- 
cies between the SLH-retrieved and sounding-based profiles for Q> are larger than 
those for heating, key features of the vertical profiles agree well. The Q 2 drying 
magnitudes are smaller than the sounding-derived magnitudes in the upper tro- 
posphere. This is because the SLH algorithm is severely limited by the inherent 
sensitivity of the PR that can detect only precipitation-sized particles. 

The SLH algorithm was also applied to PR data for February 1998 (El Nino) 
and February 1999 (La Nina). The differences in the SLH-estimated Q 2 between 
February 1998 and February 1999 are consistent with those in the SLH-estimated 
heating between February 1998 and February 1999. It is also shown that the 
SLH algorithm can estimate differences of Q 2 between the western Pacific and the 
Atlantic Ocean as well as those of heating, which are consistent with the results 
from the budget study. 
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1. Introduction 


The Tropical Rainfall Measuring Mission (TRMM; Simpson et al. 1988, 1996; 
Kummerow et al. 2000), a joint Japanese-U.S. cooperative Earth Probe satellite, 
was successfully launched in 1997 to advance understanding of the global en- 
ergy and water cycle. The TRMM satellite has been in operation for more than 9 
years, providing the distribution of rainfall throughout the Tropics using microwave 
observations from the Precipitation Radar (PR) and the TRMM Microwave Imager 
(TMI). Estimating vertical profiles of latent heating released by precipitating cloud 
systems is one of the key objectives of TRMM, together with accurately measuring 
the horizontal distribution of tropical rainfall (see a review by Tao et al. 2006). 

The PR is the first space-borne precipitation radar and can provide height infor- 
mation based upon the time delay of the precipitation-backscattered return power 
(Kozu et al. 2001 ; Okamoto 2003). This allows for vertical profiles of precipi- 
tation to be obtained directly over the global Tropics (iguchi et al. 2000). The 
classification between convective and stratiform regions of mesoscale convective 
systems (MCS) became more straightforward utilizing observed precipitation pro- 
files (Awaka et al. 1998). The accuracy of this classification is very important 
for estimating latent heating because differences in diabatic heating profiles exist 
between convective and stratiform regions of MCSs (Houze 1982; Johnson and 
Young 1 983). For convective regions of MCSs, the heating profile has warming 
at all levels with a maximum at midlevels, whereas in stratiform regions there is 
a warming peak in the upper troposphere and a cooling peak at low-levels. The 
resulting MCS heating profile is positive at all levels but with a maximum value in 
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the upper troposphere (“top-heavy” profile). 

The convective-stratiform heating (CSH) algorithm has been developed, based 
on the assumption that the shape of the overall MCS heating profile is deter- 
mined by the relative amounts of convective and mesoscale heating, which are 
proportional to the relative amounts of convective and stratiform precipitation (Tao 
et al. 1993a). An appropriate selection of model-generated convective and strat- 
iform heating profiles, normalized by the convective and stratiform rainfall, from 
the lookup table is very important for the CSH algorithm (Tao et al. 2000). Schu- 
macher et al. (2004) demonstrated the horizontal variation in the heating profile 
across the tropics calculated from TRMM PR observations using a method simi- 
lar to the CSH algorithm, except they used simpler, assumed profiles. In addition, 
they input their heating profiles into an idealized climate model to determine the re- 
sponse of the large-scale circulation to the heating patterns. However, differences 
in storm height and distribution of latent heating within a convective or stratiform 
column were not considered in thse approaches. Recently, Back and Bretherton 
(2006) suggest that these factors are critically important using reanalysis datasets. 

The concept of the spectral approach originates from Austin and Houze ( 1973 ) 
and Houze (1973) in which precipitation top heights (PTHs) observed by surface- 
based radar data were utilized in estimating the vertical mass transports (propor- 
tional to latent heating) by cumulus-scale convection as a function of storm fop 
height. This spectral approach has been extended by Houze and Leary ( 1976 ), 
Leary and Houze (1980), Houze et al. (1980) and Cheng and Houze ( 1980 ). 
Takayabu (2002) applied this spectral concept to obtain a spectral expression of 
precipitation profiles to examine convective and stratiform rain characteristics as a 
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function of PTH over the equatorial area (ICTN-ICPS) observed by the TRMM PR. 
In her study, nadir data from PR2A25 version 5 (Iguchi et al. 2000) for the period 
of 1998-1999 were utilized, and convective and stratiform precipitation were sep- 
arated based on the TRMM PR version 5 2A23 convective-stratiform separation 
algorithm. Convective rain profiles show near monotonic change with cumula- 
tive frequency. Stratiform rain profiles consist of two groups. One group consists 
of shallow stratiform rain profiles, which are very weak and increase downward. 
The other group consists of anvil rain profiles, characterized by maximum intensity 
around the melting level, much less intensity above, and a downward decrease 
below as indicated in traditional radar observations (e.g., Leary and Houze 1979). 
Schumacher and Houze (2003b) suggested that since the shallow, isolated echoes 
represent warm rain processes, they should be classified as convective. After 
this suggestion of Schumacher and Houze, the spectral plots of Takayabu (2002: 
Fig. 1) were revised by reclassifying shallow, isolated rain (rain type 15 in product 
2A-23) as convective (Shige et al. 2004: Fig. 2). 

Based on the results of the spectral precipitation statistics of Takayabu (2002), 
the Spectral Latent Heating (SLH) algorithm was developed for the TRMM PR in 
part ! of this study (Shige et al. 2004, hereafter Part I). The method uses PR in- 
formation (i.e., PTH, precipitation rates at the surface and melting level, and rain 
type) to select the heating profiles in lookup tables. Heating profile lookup tables for 
the three rain types-convective, shallow stratiform, and anvil rain (deep stratiform 
with PTH higher than the melting level)~were derived with numerical simulations 
of tropical cloud systems in TOGA COARE (Tropical Ocean Global Atmosphere 
Coupled Ocean-Atmosphere Response Experiment - Webster and Lukas 1992) 
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utilizing a cloud-resolving model (CRM). For convective and shallow stratiform re- 
gions, the lookup table is based on the PTH. Considering the sensitivity of the 
PR, we used a threshold of 0.3 mm h _1 to determine the PTH. Properties of the 
convective and shallow stratiform heating profiles show near-monotonic change 
with PTH, suggesting that the distribution of latent heating is a strong function of 
PTH. On the other hand, the PR cannot observe the PTH accurately enough for 
the anvil regions because of its insensitivity to the small ice-phase hydrometers 
(Heymsfield et al. 2000). Thus, for the anvil region, the lookup table refers to the 
precipitation rate at the melting level P m instead of PTH. The utilization of PTH and 
P m provides two distinct advantages for the SLH algorithm. First, the differences 
in heating profiles between the shallow convective stage and the deep convective 
stage can be realistically retrieved. Second, heating profiles in the decaying stage 
with no surface rain can also be retrieved. Preliminary applications of the SLH al- 
gorithm using TRMM PR data have been done. Tao et al. (2006) presented latent 
heating structure for a tropical Pacific Typhoon and a tropical, oceanic MCS esti- 
mated by the SLH algorithm. Morita et al. (2006) examined latent heating structure 
of the Madden-Julian Oscillation (MJO: Madden and Julian 1994) using the SLH 
estimates. On the other hand, Grecu and Olson (2006) used a procedure similar 
to the SLH algorithm to assign a heating profile physically consistent with each 
precipitation profile derived from the the combined TRMM PR-TMI algorithm. 

In the subsequent paper (Shige et al. 2007, hereafter Part II), the universal- 
ity of the lookup table produced from COARE simulations used in the SLH algo- 
rithm was examined for its global application to TRMM PR data and improvements 
were made to the SLH algorithm. In the revised procedure for convective heating 
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retrieval, the upper-level heating amplitude due to ice processes and lower-level 
heating amplitude due to liquid water processes are determined separately. In the 
stratiform region, the heating profile is shifted up or down by matching the melting 
level in the COARE lookup table with the observed one. The revised SLH algorithm 
was applied to PR data, and the results compared to Q x profiles derived diagnos- 
tically from SCSMEX (the 1998 South China Sea Monsoon Experiment) sounding 
data (Johnson and Ciesielski 2002), where Q x is the apparent heat source budget 
defined in Yanai et al. (1973). Key features of the vertical profiles agree well, par- 
ticularly the level of maximum heating. Differences of heating profiles between the 
western Pacific (more “top-heavy”) and the Atlantic Ocean (more “bottom-heavy”) 
estimated by the SLH algorithm are also consistent with those from the budget 
study (Thompson et al. 1979, see a review by Cotton and Anthes 1989). By 
utilizing the information about precipitation profiles, the SLH algorithm retrieves 
differences in the shape of convective heating profile between the eastern Pacific 
and the western Pacific during the cold phase. The differences in the shape of 
convective heating profiles across the Pacific are consistent with the results from 
reanalysis datasets (Back and Bretherton 2006). 

Large-scale models (i.e., general circulation an climate models) require not only 
the vertical distribution of Q x but also that of Q 2 , where Q 2 is the apparent moisture 
sink budget defined in Yanai et al. (1973). Rajendran et al. (2004) developed a 
new empirical cumulus parameterization scheme (ECPS), based on a procedure 
to improve the vertical distribution of heating and moistening over the tropics. The 
ECPS utilizes observed profiles of Q x derived using TMI-CSH algorithm (Tao etal. 
1 993a) and the European Centre for Medium-Range Weather Forecasts (ECMWF) 
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analysis based Q 2 . Comparisons between short range ECPS forecasts and those 
with the modified Kuo scheme (Krishnamurti et al. 1983), show a very marked im- 
provement in the skill in ECPS. More improvements may be achieved with ECPS 
utilizing profiles both Q x and Q 2 derived from TRMM data (Rajendran 2005, per- 
sonal communication). However, none of TRMM heating algorithms estimates Q 2 
profiles. 

In this study, the SLH algorithm is used to estimate Q 2 profile. First, differ- 
ences between Q x and Q 2 are examined. Second, Q 2 profiles are reconstructed 
from CRM-simulated parameters (i.e., PTH, precipitation rates at the surface and 
melting level, and rain type) and compared to the “true” CRM-simulated Q 2 pro- 
files, which are computed directly by the model water vapor equation. Third, the 
SLH algorithm is applied to PR data and the results compared to Q 2 profiles de- 
rived diagnostically from SCSMEX sounding data (Johnson and Ciesielski 2002 ). 
It is also applied to PR data for February 1998 and February 1999 to highlight the 
differences between El Nino-La Nina events. 


2. Approach 

Figure 1 shows the procedure for refining and validating the SLH algorithm. 
Due to the scarcity of reliable validation data and difficulties associated with the 
collocation of validation data and satellite measurements, a consistency check of 
the SLH algorithm is performed using CRM-simulated precipitation profiles as a 
proxy for the PR data. The SLH algorithm estimates Q 2 mainly due to precipita- 
tion processes, because it is severely limited by the inherent sensitivity of the PR 
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that can detect only precipitation-sized particles. Hereinafter, Q 2 estimated by the 
SLH algorithm are represented as Q 2p to bring attention to this point. Algorithm- 
reconstructed Q 2p profiles are derived from CRM-simulated precipitation profiles 
and compared to CRM-simulated “true” Q 2 profiles, which are computed directly 
from the model water vapor equation. The consistency check is a useful and nec- 
essary precondition for the application of the algorithm to actual TRMM PR data. 
CRM-simulated data from TOGA-COARE, GATE (Global Atmospheric Research 
Program Atlantic Tropical Experiment - Houze and Betts 1981), SCSMEX and 
KWAJEX (the 1 999 Kwajelin Atoll field experiment - Yuter et al. 2005) are used 
as a consistency check. Only precipitation over ocean is considered in the current 
investigation. In this paper, the SLH algorithm is applied to PR data and the results 
are compared to Q 2 profiles derived diagnostically from SCSMEX sounding data 
(Johnson and Ciesielski 2002). 

The CRM used in this study is the two-dimensional version of the Goddard Cu- 
mulus Ensemble (GCE) model and is primarily documented in Tao and Simpson 
( 1993 ). Recent improvements were presented in Tao (2003) and Tao et al. (2003). 
The model includes solar and infrared radiative transfer processes, and explicit 
cioud-radiation interactive processes (Tao et al. 1 996). Simulations presented in 
this study employ a parameterized Kessler-type two-category liquid water scheme 
(cloud water and rain) and a three-category ice-phase scheme (cloud ice, snow, 
and graupel) by Rutledge and Hobbs (1984). A “new saturation technique” (Tao 
et al. 2003 ), which allows the temperature to change after the water phase before 
the ice phase is treated, is used. This technique yields far less cloud water and 
significantly more cloud ice above the freezing level, which is much more in line 
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with in situ aircraft measurements of there being almost no super cooled cloud 
water in the anvil (Stith et al. 2002; Heymsfield et al. 2002). Subgrid-scale (turbu- 
lent) processes in the GCE model are parameterized using a scheme based on 
Klemp and Wilhelmson (1978) and Soong and Ogura (1980). The effects of both 
dry and moist processes on the generation of subgrid-scale kinetic energy have 
been incorporated in the model. The model domain is 1024 km in the x direction 
(horizontal) and 22.4 km in the z direction (vertical). The horizontal resolution is 
1000 m. The vertical resolution varies from 100 m at the lower boundary to 1000 
m at the top of the domain. The time step is 12 s. 

In this study, tropical convective systems in TOGA-COARE, GATE, SCSMEX 
and KWAJEX are simulated with an approach, so-called cloud ensemble model- 
ing. In this approach, many clouds of different sizes in various stages of their 
lifecycles can be present at any model simulation time. Observed large-scale ad- 
vective tendencies of temperature, moisture, and horizontal momentum are used 
as the main large-scale forcings that govern the GCE model in a semi-prognostic 
manner (Soong and Ogura 1980). These are applied uniformly over the model 
domain with the assumption that the model domain is considerably smaller than 
the large-scale disturbance. Large-scale advective tendencies for temperature T 
and specific humidity q are defined as: 



V 0 bs ' T" i'abs ^obs 


dT obs 

dp 


+ 


®obs 

^obs 

Op 


“V obs 


^ Qobs W obs 


dq 0 bs 

dp 


( 1 ) 

(2) 


and were derived from sounding networks deployed during TOGA COARE, GATE, 
SCSMEX and KWAJEX. Here v is the horizontal wind vector, u the vertical pres- 


8 



sure velocity, a the specific volume, and C v the heat capacity at constant pressure. 

Since accurate calculations of the large-scale horizontal momentum forcing 
terms are difficult to obtain from observations in the Tropics (Soong and Tao 1 984), 
the terms are instead replaced by a nudging term: 


(dv\ V-Vobs 

- 


(3) 


where v is the model domain averaged horizontal velocity, v ob3 the observed large- 
scale horizontal vector over the sounding networks, and r the specified adjustment 
time scale of six hours. This method constrains the domain-averaged horizontal 
velocities to follow the observed values, and thereby provides a simple means in 
controlling the cloud system dynamics by the large-scale momentum and shear. 
Cyclic lateral boundary conditions are incorporated to ensure that there is no ad- 
ditional heat and moisture forcing inside the domain other than the imposed large- 
scale forcing. • 

The accuracy of the convective-stratiform separation affects the inference of 
the vertical distribution of heating and moistening. The TRMM PR rain-type clas- 
sifications, in which bright band identification is very important, cannot be directly 
applied to GCE outputs (Awaka et al. 1996). The microphysical schemes utilized 
in CRMs (e.g., Lin et al. 1983; Rutledge and Hobbs 1984) typically do not contain 
an explicit description of the partially melted precipitation particles that lead to a 
bright band of enhanced radar reflectivity. Thus, the GCE convective and strati- 
form separation method (Lang et al. 2003) is used with some modifications done 
in Part I to maintain the consistency with the TRMM PR rain-type classification. 
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3. Heat and moisture budgets 


In diagnostic studies (Yanai et al. 1973; Yanai and Johnson 1993), it is custom- 
ary to define the apparent heat source Q x and the apparent moisture sink Q 2 of a 
large-scale system by averaging horizontally the thermodynamic and water vapor 
equations as 


Q 1 = 7T 


'do 

dt 


+ V • V9 + w 


dO' 

dz 


(4) 


Q2 — ~ 


Ly dQy 

cr p [m 


+ V • V q v + w 


dq v 

dz 


(5) 


where n = (p/P 0 o) R/Cp is the nondimensional pressure, P 00 the reference pressure 
(1000 mb), C p the specific heat of dry air at constant pressure, and R the gas 
constant for dry air. 

Q 1 and Q 2 can be directly related to the contributions of cloud effects, which 
can be explicitly estimated by CRMs as: 



1 dpw'O' 
p dz 


— v' • VO' 4 - Dq 


+ LH + Qr 1 


Q2 


L v 

a 


1 dpw’q' v 
p dz 


v' • Vq' v + D, 


C[v 


+ NC. 


( 6 ) 

(7) 


The overbars denote horizontal averages; the primes indicate deviations from the 
horizontal averages, p is the air density, and Q R the cooling/heating rate associated 
with radiative processes. The subgrid-scale (smaller than the cloud scale) diffusion 
are represented by D e and D gv that are usually small compared to other terms 
above the boundary layer (Soong and Tao 1980). The term LH and NC are the 
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net latent heating and the net condensation due to the phase change of water: 

+ + ( 8 ) 

Op L/p Op 

NC= ~(c~e + d~s). (9) 

Op 

L v , L f and L s are the latent heats of vaporization, fusion and sublimation, re- 
spectively. Variables c, e, d and 5 stand for the rates of condensation of cloud 
droplets, evaporation of cloud droplets and raindrops, freezing of cloud droplets 
and raindrops, melting of ice crystals, snow and graupel, deposition of ice crys- 
tals and sublimation of all ice hydrometeors, respectively. These processes are 
not directly detectable with remote sensing (or for that matter, with in situ mea- 
surements). Thus, heating and moistening retrieval schemes depend heavily on 
the use of CRM. The first terms on the right-hand side of Eqs (6) and (7) are the 
vertical eddy heat and moisture flux convergences from upward and downward 
cloud-scale motions, respectively. The second term are the horizontal eddy heat 
and moisture flux convergences, respectively. 

We accumulate Qi and Q 2 over a period of 5 minutes for each data sampling, 
since accumulation over long periods are inadequate for growing convective cells 
(Shige and Satomura 2000: Fig. 4a) and moving convective systems. Additional 
sensitivity tests with periods of 1 minute and 2 minutes indicate that there are more 
quantitative differences in convective and stratiform components of Q 2 than those 
of Qi. In this study, the two-dimensional version of the GCE model was used. 
Real clouds and cloud systems are three-dimensional. Grabowski et al. (1998) 
and Donner et al. (1999) found larger temporal variability in the two-dimensional 
simulation than in the three-dimensional simulation. Thus, the sensitivity of Q 2 in 
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convective and stratiform regions to accmulated time might be smaller in the three- 
dimensional simulation than in the two-dimensional simulation, but it is beyond the 
scope of this study. Hereafter Q x and Q 2 accumulated over a period of 5 minutes 
for each data sampling are represented as instantaneous Q x and Q 2 to maintain 
consistency with Part I and II. 

Figures 2a-f show GCE-simulated average profiles of heat budget [LH, Q x -Q R 
(hereinafter Q 1R ) without and with horizontal eddy heat flux] and moisture budget 
[NC, Q 2 without and with horizontal eddy heat flux] for the TOGA COARE (19-26 
December 1992) case, respectively. The difference between LH and NC around 
4.4 km is large (Figs 2a, d). The maximum of NC around 4.4 km can be explained 
by a decrease of saturation water vapor mixing ratio, which is due to the melting 
cooling (Guichard et al. 1997). It is also evident that the NC magnitudes are 
smaller than the LH magnitude above 6 km, because freezing is included in LH 
(Eq. (8)), but not in NC (Eq. (9)). 

The profile of Q 1R for the total region (Fig. 2b) is close to that of LH for the 
total region (Fig. 2a) except the melting level where the vertical eddy heat flux 
convergence compensates for the distinct LH cooling due to the melting for the 
total region (Sui et al. 1994; Shie et al. 2003). On the other hand, the profile of Q 2 
for the total region (Fig. 2e) does not follow that of NC for the total region (Fig. 2d) 
and has a more complicated shape. The vertical eddy moisture flux convergence is 
a major contributor to the Q 2 while the contribution to the Q XR budget by the vertical 
eddy heat flux convergence is minor, as noted by previous studies (Soong and Tao 
1980; Tao and Soong 1986; Lafore et al. 1988; Chong and Hauser 1990; Simpson 
and Tao 1993; Tao et al. 1993b; Caniaux et al. 1994; Guichard et al. 1997). This 
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difference is explained by the fact that the vertical gradient of the mean water vapor 
is larger than that of the mean potential temperature. 

The horizontal eddy heat and moisture flux convergences have been neglected 
in most of the previous studies. Neglecting the horizontal eddy heat and moisture 
flux convergences is appropriate for Q 1R and Q 2 for the total region because differ- 
ences between profiles of Qw and Q 2 without and with horizontal eddy heat flux 
are negligible. However, as discussed in Part II, there are the differences between 
profiles of Q XR without and with horizontal eddy heat flux in the convective and 
stratiform regions (Figs 2b, c), and thus horizontal eddy heat flux should not be 
neglected. Larger differences between profiles of Q 2 without and with horizontal 
eddy moisture flux in the convective and stratiform regions can be found (Figs 2e, 
f). Larger differences of Q 2 than Q lR can be attributed to the greater horizontal 
variability of moisture than temperature. 

Lafore et al. (1988) and Caniaux et al. (1994) argued that the double-peak 
structure in Q 2 often seen in tropical budget studies (Reed and Recker 1971 ; Yanai 
et al. 1973; Johnson 1976) is a result of midtropospheric moistening by vertical 
eddy moisture convergence in the convective region, while Johnson (1984) ex- 
plained this structure as a consequence of the combined but vertically separated 
drying effects of convective and stratiform regions. Although the TOGA-COARE 
I FA (Intensive Flux Array) Q 2 profile does not exhibit a prominent double-peak 
structure as shown in Figs 2f, there are times when the TOGA-COARE IFA Q 2 pro- 
files do resemble those obtained by the previous study (Johnson and Lin 1997). 
An example is during a period of heavy rainfall on 11-17 December (Fig. 3; see 
also Fig. 3 of Johnson and Lin 1997). A double peak in Q 2 is present (though not 
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very pronounced). Profiles of Q 2 without horizontal eddy heat flux in the convective 
and stratiform regions (Figs 2b) are in good agreement with those of Q 2 in the con- 
vective and stratiform regions simulated by Caniaux et al. (1994: their Fig. 22b).. 
Caniaux et al. (1994) neglected horizontal eddy moisture flux convergence, be- 
cause it is one order of magnitude smaller than NC and vertical eddy moisture 
flux convergence. This is true for the total region, but not for the convective and 
stratiform regions. The horizontal eddy moisture flux convergence compensates 
midtropospheric moistening by vertical eddy moisture convergence in the convec- 
tive regions during this period (Fig. 3c), as well as the period of heavy rainfall on 
19-26 December (Fig. 2f). From our results, the double-peak structure is likely ex- 
plained as a consequence of the combined but vertically separated drying effects 
of convective and stratiform regions. It should be noted that Lafore et al. (1988) 
and Caniaux et al. (1994) examined a tropical continental squall line. Updrafts 
associated with tropical continental cases are expected to be stronger than those 
associated with tropical oceanic ones. Thus, even if the horizontal eddy moisture 
flux convergence is included, their results might be different from ours. 

The SLH algorithm estimates Q lR and Q 2 mainly due to precipitation processes 
{Quip and Q 2p ), because it is severely limited by the inherent sensitivity of the PR 
that can detect only precipitation-sized particles. However, it is noticed in Figs 2 
and 3 that the Q 2 in the nonprecipitating region is the same order as that in the 
stratiform region, while the Q iR profile in the nonprecipitating region is negligible. 
Thus, we cannot expect that differences between Q 2p and Q 2 are as small as those 
between Q XRp and Q 1R . 


14 



4. Algorithm 


a. Construction of lookup tables 


The Q 2p profile lookup tables have been constructed for the three rain types; 
convective, shallow stratiform, and deep stratiform with and without surface rain 
(Fig. 4). For construction of lookup tables, the GCE-simulated precipitation profiles 
and corresponding Q 2p profiles from the four sub-periods of 9-day durations (10-18 
December 1992, 27 December 1992 -4 January 1993, 9-17 February 1993, and 
18-26 February 1993) are used. 

Figure 4a shows the lookup table for convective rain. The GCE-simulated pre- 
cipitation profiles with 0.3 mm h _1 precipitation top threshold and corresponding 
Q 2p profiles are accumulated and averaged for each PTH with model grid inter- 
vals. Properties of convective Q 2p profiles show less monotonic changes with PTH 
than Qi Rp do (Part II; Fig. 5a) because contribution of eddy fluxes is larger in Q 2p 
than in Q lRp . 

As for the stratiform precipitation, the GCE-simulated precipitation profiles with 
0.3 mm h _1 precipitation top threshold and corresponding Q 2p profiles are also 
averaged for each PTH (Fig. 4b). However, as mentioned earlier, we cannot ex- 
pect the PR to retrieve the small ice-phase precipitation, and observe the PTH 
accurately enough in the upper-level regions of the anvils where small ice-phase 
hydrometers dominate as mentioned earlier. Therefore, only the stratiform Q 2p pro- 
files with PTH lower than 4.4 km (shown in Fig. 4b), characterized by moistening 
around PTH levels, are used as the shallow stratiform heating profile lookup table 


Fig. 4 
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for the SLH algorithm. 

Figure 4c shows lookup table for anvil (deep stratiform with PTH higher than 
the melting level) rain. PR can measure the precipitation rate at the melting level 
as well as surface-based radar ( e.g . Leary and Houze 1979), although it cannot 
observe the PTH accurately enough in the upper-level regions of the anvils. Thus, 
for the anvil region, the lookup table refers to the precipitation rate at the melting 
level P m instead of PTH (Part I). The anvil profiles with PTH higher than the melting 
level are characterized by upper-level drying and lower-level moistening, which is 
also found in observations {e.g. Johnson and Young 1983). The upper ievel drying 
in these anvil regions is largely due to condensation and deposition, while the 
lower-level moistening is largely due to evaporation of raindrops. 


b. Procedure of Q 2p retrieval 


Procedure of Q 2p retrieval is the same as that of heating retrieval (Part I, II) ex- 
cept for using the Q 2p tables (Fig. 4). The SLH algorithm utilizes the observed infor- 
mation on precipitation type (convective or stratiform), PTH, P s , and P m . The con- 
vective region is then separated into a shallower convective region and a deeper 
convective region. On the other hand, the stratiform region is separated into a 
shallow stratiform region and an anvil region based on the PTH in comparison with 
the melting level. The algorithm then derives Q p 2 profiles based on the following 
procedure for each of the four regions (Fig. 5). 

For shallower convective regions as well as shallow stratiform regions, a Q 2p 
profile corresponding to the PTH is selected in the convective Q 2p profile (Fig. 4a) 


Fig. 5 


16 




and PTH-based shallow stratiform Q 2p profile (Fig. 4b) lookup tables, respectively. 
The amplitude is determined by 

Q(z) = %±-P.. (10) 

* s 

For deeper convective regions, the upper-level drying due to ice processes is 

determined by 

= (id 

where P f is the precipitation rate at the level separating upper-level drying from 
lower-level drying. Likewise, the lower-level drying due to liquid water processes 

is determined by 

Q(z) lm , = • P.. (12) 

For anvil regions with PTH higher than the melting level, on the other hand, 
the Q 2p profile corresponding to P m is selected in the anvil Q 2p profile lookup table 
(Fig. 4c). The upper-level drying amplitude is then determined by 

Qm 9 h{*) = . p m , (13) 

m 

Evaporative moistening rate below the melting level in the anvil regions is propor- 
tional to the reduction of precipitation profile toward the surface from the melting 
level. The downward decrease of the intensity of anvil rain below the melting layer 
has been statistically indicated in the PR observation (see Fig. 2 of Takayabu, 
2002) as well as in case studies of traditional radar observations (e.g. Leary and 
Houze 1979). Thus, the algorithm computes the lower-level moistening amplitude 
Q to . w as a function of the difference of rainfall rate between the surface and the 
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melting level P m — P s : 


Qloviiz) = • (Pm - Ps) (14) 

4 m -L s 

The above procedure allows us to estimate the Q 2p profiles in the anvil regions 
both with and without surface precipitation. The Q 2p profile is shifted up or down 
by matching the melting level of the COARE lookup table with the observed one. 


5. Consistency check 

a. Time averaged profiles 


Four eight-day periods from TOGA-COARE (19-26 December 1992), GATE 
(Sep 1-8 1974), SCSMEX (Jun 2-9 1998) and KWAJEX (Sep 6-13 1999) are used 
for a consistency check of the SLH algorithm, as shown in Fig. 6. For each period, 
Q 2p profiles were reconstructed using the simulated parameters (i.e. PTH, convec- 
tive/stratiform characteristics, P s , and P m ) as input. The algorithm-reconstructed 
Q 2 P profiles from the GCE-simulated precipitation profiles are compared with GCE- 
simulated true Q 2 profiles for the convective, stratiform, non-precipitating and total 
regions. 

For the COARE case, the SLH algorithm with the COARE lookup table pro- 
duces good agreement between the SLH-reconstructed Q 2p and GCE-simulated 
Q 2 profiles for convective and stratiform regions (Fig. 6a). The reconstructed to- 
tal Q 2p profile is in good agreement with the model except for lowest-levels where 
moistening due to non-precipitating processes is dominant. 


Fig. 6 
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For the GATE case, the COARE lookup table results in less agreement between 
the SLH-reconstructed and GCE-simulated profiles for convective and stratiform 
regions (Fig. 6b). The SLH-reconstructed drying at z = 1-3 km is weaker than the 
GCE-simulated for the convective region, while the SLH-reconstructed moistening 
at z = 1-3 km is weaker than the GCE-simulated for the stratiform region. In 
spite of compensating errors at z = 1-3 km from each component (convective and 
stratiform), the SLH-reconstructed drying at z = 2-6 km is stronger than the GCE- 
simulated for the total region. This discrepancy is mostly explained by moistening 
for non-precipitating region that SLH cannot reconstruct. 

The COARE lookup table produces good agreement between reconstructed 
Q 2p and simulated Q 2 profiles for the SCSMEX convective and stratiform regions 
(Fig. 6c). In addition to good agreement for the convective and stratiform regions, 
contribution of non-precipitating region to the SCSMEX total Q 2 profile is small, 
leading good agreement reconstructed Q 2p and simulated Q 2 profiles for the SC- 
SMEX total region. 

For the KWAJEX case, the SLH-reconstructed drying below z = 4 km is weaker 
than the GCE-simulated for the convective Q 2 profiles, and while the SLH-reconstructed 
moistening below z = 4 km is weaker than the GCE-simulated for the stratiform Q 2 
profiles, similar to the GATE case. Because of compensating errors below z = 4 km 
from each component (convective and stratiform), the reconstructed total Q 2 pro- 
file is in good agreement with the model except for lowest-levels where moistening 
due to non-precipitating processes is dominant. 

The differences between the SLH-reconstructed Q 2p and GCE-simulated Q 2 
shown here are larger than those between the SLH-reconstructed Qm p and GCE- 
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simulated Q 1R (Part II: Fig. 9). This is because the Q 2 profiles are noisier than 
those of Q lR due to much more contribution of eddy fluxes. The SLH algorithm 
cannot estimate the effect (moistening) of non-precipitating area, leading to a dis- 
crepancy with the GCE results. Nevertheless, the SLH-reconstructed total Q 2p pro- 
files are in good agreement with the GCE-simulated ones, espcially for the case 
where contribution of non-precipitating region to the Q 2 budget is small. 

b. Error estimation 

In the SLH algorithm, lookup tables are constructed based on the assumption 
that heating profiles correspond statistically to precipitation profiles or precipitation 
parameters (i.e., PTH, P m ). However, the instantaneous grid cell relationship be- 
tween precipitation profiles and Q 2p is somewhat ambiguous. Part II performed a 
preliminary evaluation of the horizontally averaged estimates and found that hori- 
zontal averaging over —30 km width was required to reduce random errors in the 
SLH-reconstructed heating profiles to acceptable levels. 

Following Part II, a preliminary evaluation of the horizontally-averaged esti- 
mates of Q 2p for the COARE, GATE, SCSMEX and KWAJEX periods used in the 
consistency check is performed. Q 2p profiles were reconstructed grid by grid for 
each 8-day periods using the simulated parameters as input. Then, the differences 
between the reconstructed Q p2 profiles and the simulated Q 2 were examined sta- 
tistically to see the errors in the instantaneous grid cell estimates using the table 
method. Larger root-mean-square (rms) errors were found for Q 2p than with Q lRp 
(Part II: Fig. 10). While the rms errors of Q lRp at the PR foot-print scale (4 km) 
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are smaller than 16 K 2 h~ 2 , those of Q 2p are much larger than 16 K 2 h 2 . Larger 
contribution of eddy fluxes in the Q 2p budget than in the Q\r p budget is a major 
reason. Averaging over ~60 km in width reduces the rms to about 1 K 2 h~ 2 . From 
these results, averaging over ~60 km in width is recommended in order to use the 
SLH algorithm estimates of Q 2 quantitatively. 


6. PR applications 



In this section, the revised SLH algorithm is applied to precipitation profiles 
from version 6 of the TRMM PR 2A25 data set, which is instantaneous and at 
footprint-scale (i.e., a level-2 product). Recently, Shige et al. (2006) investigated 
the consistency between TMI-observed brightness temperatures (TBs) and those 
simulated from PR2A25 V5 and V6 rain profiles using a radiative transfer model. 
They showed that simulated TBs from PR V6 exhibits better agreement with ob- 
served ones than those from PR V5, implying the algorithm improvements. Bright- 
band height estimates from version 6 of TRMM PR 3A25, gridded five degree spa- 
tial resolution monthly composite of instantaneous and footprint-scale data (PR 
2A25), are also used as the melting levels. 


a. Comparison of Q-> profiles over the SCSMEX NESA region 


in Part II, the accuracy of the SLH-retrieved heating was evaluated by compar- 
ing with a rawinsonde-based analysis of diabatic heating for the SCSMEX NESA 
derived by Johnson and Ciesielski (2002). Key features of the vertical profiles 
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agree well, particularly the level of maximum heating. Magagi and Barros ( 2004 ) 
and Grecu and Olson (2006) also compared their results against heating estimates 
over the SCSMEX NESA. Here, the accuracy of the SLH-retrieved Q 2p are evalu- 
ated by comparing with a rawinsonde-based analysis of Q 2 for the SCSMEX NESA 
derived by Johnson and Ciesielski (2002). 

Figure 8 shows a comparison between SLH-retrieved Q 2p from version 6 of 
the TRMM PR data sets and sounding-based Q 2 during the campaign’s most con- 
vectively active period (May 15 - Jun 20 1998). Mapes et al. (2003) suggested 
that averages of about 30 days reduce sampling errors in the rainfall rate estimate 
(proportional to integrated Q x or Q 2 ) to 1 0% for the SCSMEX NESA. There is good 
agreement in several key features of the vertical profiles. It is evident from Fig. 8 
that the Q 2p drying magnitudes are smaller than the sounding-derived magnitudes 
in the upper troposphere. The Q lRp heating magnitudes were also smaller than the 
sounding-derived magnitudes in the upper troposphere (Part II: Fig. 11). The SLH 
algorithm is severely limited by the inherent sensitivity of the PR that can detect 
only precipitation-sized particles. During the growing phase of a congestus cloud, 
cloud top and radar-echo top may correlate well (Kingsmill and Wakimoto 1991). 
However, during the decaying phase of a cumulonimbus cloud and in stratiform 
regions, the two tops may differ significantly, leading to underestimate of drying in 
the upper troposphere. Therefore, measurements from other sensors [e.g., VIRS 
(Visible and Infrared Scanner)] will have to be integrated to obtain a more accurate 
estimation of moistening profiles, but it is beyond the scope of this study. 
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b. El Nino-La Nina 


Figures 9a, b show the monthly mean surface rainfall (mm day -1 ) for February 
1998 and February 1999, respectively. Drying/moistening structures over the six 
oceanic regions (western Pacific, central Pacific, east Pacific, south Pacific, In- 
dian Ocean and Atlantic Ocean) shown in Fig. 9 will be examined and compared. 
Heating structures over the same regions were examined in Part II. 

Figures 10 and 1 1 show the monthly mean convective, stratiform and total Q 2p 
profiles derived from the SLH algorithm for six locations over the tropical oceans 
for February of 1998 and February of 1999, respectively. Also total Q XRp profiles 
derived from the SLH algorithm are shown. 

A maximum at low-levels (~0.5 km) is found in the SLH-estimated total Q 2p pro- 
file over the western Pacific for February 1998 (Fig. 10a). This low-level maximum 
in the SLH-estimated total Q 2p profile comes from the SLH-estimated convective 
Q-) P profile with a low-level maximum, reflecting the abundance of shallow convec- 
tion. Diagnostic budget studies over west Pacific regions (Reed and Recker 1971 ; 
Nitta 1972; Yanai et al. 1973; Lin and Johnson 1996) indicate a double-peak struc- 
ture with a minimum near 4 km. While the SLH-estimated mean Q 2p profile for 
February 1999 show a minimum near 4 km (tough not very pronounced), the SLH- 
estimated mean Q 2p profile with the low-level maximum for February 1998 does 
not resemble those determined from the diagnostic budget studies. It should be 
noted that diagnostic budget studies over the western Pacific do not contain peri- 
ods corresponding to the warm phases of ENSO, except for two months out of the 
period from March to July of 1958 in Nitta (1972). The SLH-estimated mean Q 2p 


Fig. 9 


Fig. 10 
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profile for February 1 998 resembles the mean Q 2 profile with a lower-level maxi- 
mum over the GATE region (Nitta 1 978; Thompson et al. 1 979), which is due to the 
different cloud population and lower SSTs in the GATE region. Deep convection 
over the western Pacific is suppressed during the warm phase of ENSO (February 
1998) relative to the cold phase (February 1999) due to lower sea surface tem- 
peratures. Thus, the difference in the SLH-estimated mean Q 2p profile between 
February 1998 and February 1999 may be reasonable. 

There are dramatic differences in total Q 2p profiles over the central and eastern 
Pacific between the 1998 El Nino event and the 1999 La Nina event. First, these 
differences can be attributed to those in stratiform rain fraction (Schumacher and 
Houze 2003a). Second, the shape of convective Q 2p profiles also affect total Q 2p 
profiles. The convective Q 2p profiles are shallower during La Nina than El Nino. 

Part II pointed out that the SLH-estimated heating profile over the south Pa- 
cific for February 1998 is very similar to the vertical distribution of heating during 
the undisturbed BOMEX (the Barbados Oceanographic and Meteorological Exper- 
iment) period in the trade wind belts (Nitta and Esbensen 1974) and that during 
episodic trade wind regimes over the western Pacific (Johnson and Lin 1997). On 
the other hand, the SLH-estimated Q 2p profile with low-level drying over the south 
Pacific for February 1 998 is not similar to the vertical distribution of Q 2 with low- 
level moistening during the undisturbed BOMEX and that during episodic trade 
wind regimes over the western Pacific. Possible reason for the differences is that 
contribution of non-precipitating region to the Q 2 budget which SLH cannot retrieve 
is large. Another reason is that there is larger spatial variability in Q 2 profiles than 
in Qi profiles. For example, Lin and Johnson (1996) indicated that although the 
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heating profiles over different regions of the western Pacific warm Pool are similar 
to one another, the drying/moistening profiles are significantly different. 

The SLHI-estimated Q 2p profiles over the Atlantic Ocean for February 1 998 and 
February 1999 resemble the mean Q 2 profile with a lower-level maximum that was 
determined from a diagnostic budget study during GATE (Nitta 1978; Thompson 
et al. 1979) and simulated by the GCE model (see Fig. 6b). Thus, it is shown that 
the SLH algorithm can estimate differences of Q 2p between the western Pacific 
and the Atlantic Ocean as well as those of Qi, which are also consistent with the 
results from the budget study (Thompson et al. 1979, see a review by Cotton and 
Anthes 1989). 

Over the Indian Ocean, even indirect validation using the results of diagnostic 
studies for different years is difficult because of the lack of observations*. The SLH- 
estimated Q 2p profile shape over the Indian Ocean for February 1998 resemble 
that over the western Pacific for February 1999, while that over the Indian Ocean 
for February 1999 resemble that over the Atlantic Ocean. The differences in the 
SLH-estimated Q 2p over the Indian Ocean between February 1998 and February 
1999 are consistent with those in the SLH-estimated Q iR over the Indian Ocean 
between February 1998 and February 1999. 


"•Nitta (1980) did preliminary budget computations over the bay of Bengal for a very short time 

(three-day period). 
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c. Variability of Q 2p profile 


The average profiles shown in Figs. 10 and 1 1 may not be very representative, 
because there are large spatial and temporal variabilities in rainfall over a larger 
region such as those selected in Fig. 9. Following Part II, we used contoured- 
frequency-by-altitude diagrams (CFADs; Yuter and Houze 1 995) in order to provide 
actual variability of drying/moistening profile shapes. 

Many of the negative values in the low to middle troposphere in Fig. I2a-e and 
Fig. 13a, b and d are associated with evaporation in the stratiform region, suggest- 
ing large variability in stratiform rain fraction over the regions. The distinct peak in 
the frequency of drying at levels below 2 km is also evident for all six geographic 
areas for both February 1998 and February 1999 in the CFADs (Figs. 12 and 13), 
corresponding to a drying peak at lowest-levels seen in the mean convective heat- 
ing profiles (Figs. 10 and 11). It is inferred from the convective Q^-profile lookup 
table (Fig. 4a) that shallow convection with PTHs lower than 4 km accounts for this 
distinct peak. 

The CFADs of IFA-mean Q x and Q 2 during TOGA COARE have been shown 
by Johnson and Ciesielski (2000: Fig. 3). Direct comparison of the CFADs shown 
in Figs 12 and 13 against those in Johnson and Ciesielski (2000) cannot be done, 
because the CFADs shown in Figs 12 and 13 are derived form monthly mean 
Q 2p profiles at 0.5° resolution. Still, there is a striking difference. In the CFAD 
of Q 2 in Johnson and Ciesielski (2000), there is some contribution to negative 
values of Q 2 in the lowest 2 km from shallow cumulus clouds and subcloud-layer 
eddies during non-precipitating periods of the TOGA-COARE IOP (Johnson and 


Fig. 12 
Fig. 13 
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Lin 1997). On the other hand, there is almost no contribution to negative values 
of Q 2p in the lowest 2 km from shallow cumulus clouds and subcloud-layer eddies, 
although there is some contribution to negative values of Q 2p in the lowest 2 km 
from stratiform rain. Again, the SLH algorithm cannot retrieve moistening for non- 
precipitating regiondue to the sensitivity of the PR. 

7 . Summary and future work 

In this study, the SLH algorithm was used to estimate Q 2 profiles from PR data 
together with Q 1R profiles. There are two differences between Q lR and Q 2 . First, 
the eddy moisture flux convergence is a major contributor to the Q 2 while the con- 
tribution to the Qir budget by the eddy heat flux convergence is minor. Second, 
the Q 2 in the nonprecipitating region is the same order as that in the stratiform 
region, while the Q 1R profile in the nonprecipitating region is negligible. These two 
factors lead to larger errors of Q 2 than Q 1R . 

The Q 2p profiles were reconstructed from CRM-simulated parameters (i.e. PTH, 
precipitation rate at the melting level, rain rate and type) with the COARE table and 
compared then to CRM-simulated “true” Q 2 profiles, which were computed directly 
from the model water vapor equation. COARE, GATE, SCSMEX, and KWAJEX pe- 
riods were used for the consistency check. The consistency check indicates that 
discrepancies between the SLH-reconstructed Q 2p and GCE-simulated Q 2 profiles, 
especially at low levels, are larger than those between the SLH-reconstructed Q 1Rp 
and GCE-simulated Q m profiles. Larger discrepancies in Q 2p at low levels are due 
to moistening for non-precipitating region that SLH cannot reconstruct. 
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The SLH algorithm was applied to PR data, and the results compared to Q 2 
profiles derived diagnostically from SCSMEX sounding data. Although discrep- 
ancies between the SLH-retrieved and sounding-based profiles for Q 2 are larger 
than those for Q m , key features of the vertical profiles agree well. The Q 2p drying 
magnitudes are smaller than the sounding-derived magnitudes in the upper tro- 
posphere. This is because the SLH algorithm is severely limited by the inherent 
sensitivity of the PR that can detect only precipitation-sized particles. During the 
decaying phase of a cumulonimbus cloud and in stratiform regions, cloud top and 
radar-echo top may differ significantly, leading to underestimate of drying in the 
upper troposphere. 

The SLH algorithm was also applied to PR data for February 1998 (El Nino) 
and February 1999 (La Nina). The differences in the SLH-estimated Q 2 between 
February 1 998 and February 1 999 are consistent with those in the SLH-estimated 
Qir between February 1998 and February 1999. It is shown that the SLH algo- 
rithm can estimate differences of Q 2p between the western Pacific and the Atlantic 
Ocean as well as those of Q lp , which are consistent with the results from the bud- 
get study (Thompson et al. 1979, see a review by Cotton and Anthes 1989). 

Although the SLH-estimated heating profile over the south Pacific for February 
1998 is very similar to the vertical distribution of heating during the undisturbed 
BOMEX period in the trade wind belts (Nitta and Esbensen 1974) and that during 
episodic trade wind regimes over the western Pacific (Johnson and Lin 1997), the 
SLH-estimated Q 2p profile with low-level drying over the south Pacific for February 
1998 is not similar to the vertical distribution of Q 2 with low-level moistening during 
the undisturbed BOMEX and that during episodic trade wind regimes over the 
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western Pacific. Possible reason is that SLH cannot retrieve moistening for non- 
precipitating region due to the sensitivity of the PR. Measurements from other 
sensors (e.g., VIRS) will have to be integrated to obtain a more accurate estimation 

of moistening profiles. 

In this study, the two-dimensional version of the GCE model was used. Real 
clouds and cloud systems are three-dimensional. The availability of exponen- 
tially increasing computer capabilities has resulted in three-dimensional CRM sim- 
ulations for multiday periods with large horizontal domains becoming increasing 
prevalent {e.g., Grabowski et al. 1998, Donner et al. 1999). Grabowski et al. 
( 1998 ) showed the similarity in statistical aspects between the two- and three- 
dimensional CRM simulations. The reason for the similarity between the two- and 
three-dimensional CRM simulations is that the same observed large-scale advec- 
tive tendencies temperature and water vapor mixing ratio were used as the main 
forcing (Tao et al. 1987). However, there are notable differences between the 
two- and three-dimensional CRM simulations. Larger temporal variability in the 
two-dimensional simulation than in the three-dimensional simulation was found 
in Grabowski et al. (1998) and Donner et al. (1999). Tao et al. (2000, 2003) 
pointed out that the GCE three-dimensional model simulated water vapor (Q 2 ) 
budget is in better agreement with observations in the lower-troposphere than its 
two-dimensional counterpart. We are now performing the three-dimensional CRM 
simulations in order to compare look-up tables from the two- and three-dimensional 
CRM simulations. 
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Figure captions 


Fig.i. Diagram showing the procedure for refining and validating the 
spectral latent heating (SLH) algorithm. Letters denote convec- 
tive and stratiform classification (C/S), precipitation top height 
(PTH), the precipitation rate at the lowest observable level (P s ), 
and the precipitation rate at the melting level (P m ), The “?” 
means to compare and examine the Q 2 profiles reconstructed 
from the SLH algorithm with the heating profiles from modes 
simulations. This consistency check is a necessary precon- 
dition for the application of the algorithm to actual TRMM PR 
data. 

Fig. 2. Eight-day average profiles of the GCE-simulated heat budget 
[(a) LH, and Q 1R (b) without and (c) with horizontal eddy heat 
flux] and moisture budget [(d) NC , and Q 2 (b) without and (c) 
with horizontal eddy heat flux] for the total (solid), convective 
(dashed), stratiform (dotted) and non-precipitating (thin solid) re- 
gions for the TOGA COARE (19-26 December 1992) case. 

Fig. 3. Eight-day average profiles of the GCE-simulated moisture bud- 
get [(a) NC, and Q 2 (b) without and (c) with horizontal eddy heat 
flux] for the total (solid), convective (dashed), stratiform (dotted) 
and non-precipitating (thin solid) regions for the TOGA COARE 
(11-17 December 1992) case. 

Fig. 4. Ensemble-mean, GCE-simulated Q 2p profiles, plotted as func- 
tions of precipitation top height (PTH) from convective (a) and 
stratiform (b) regions, and precipitation rates at the melting level 



from anvil regions (c). Contours indicate values of confidence 
interval for the mean at the 95 % level with Student’s-t test. 
Contour interval is 2.0 K h -1 for convective and 1.0 K hr 1 for 
stratiform and anvil. Thresholds of 0.3 mm fr 1 are used for the 
precipitation top detection. 

Fig. 5 . Diagram showing the procedure for deriving latent heating pro- 
files using the spectral latent heating (SLH) algorithm. See the 
text for details. 

Fig. 6. Eight-day averaged profiles of Q 2p reconstructed by the SLH 
algorithm with the COARE lookup table (thick solid line) and Q 2 
simulated by the GCE model (dotted line) for the (a) COARE 
( 19-26 Dec 1992) case, (b) GATE (Sep 1-8 1974) case, (c) 
SCSMEX (Jun 2-9 1998) and (d) KWAJEX (Sep 6-13 1999), 
respectively. Leftmost panels are for the convective regions, 
center-left panels for the stratiform regions, center-right panels 
for the non-precipitating regions and the rightmost panels are for 
the total regions. Thin solid lines indicate differences between 
the SLH-reconstructed and the GCE-simulated profiles. 

Fig.7. The rms error in the horizontal averaged profiles between the 
SLH algorithm-reconstructed Q 2p and the GCE-simulated Q 2 
for the (a) COARE, (b) GATE, (c) SCSMEX, and (d) KWAJEX 
cases. 

Fig.8. Heating from diagnostic calculations (Johnson and Ciesielski 
2002 ) and the SLH2 algorithm using version 6 of the TRMM PR 
data sets for SCSMEX (1 5 May - 20 June 1 998). 



Fig.9. Monthly mean rainfall (mm day -1 ) derived from PR2A25 Ver- 
sion 6 for (a) February 1998 and (b) February 1999. The 
drying/moistening profiles will be compared and examined for 
the various geographic locations identified by the boxes. [From 
Shige et al. (2007)] 

Fig. io. Monthly (February 1998) mean total, convective and strati- 
form Q 2p profiles derived from the SLH2 algorithm for various 
locations. Total heating profiles derived from the CSH algorithm 
are also shown. The geographic areas are the (a) western Pa- 
cific, (b) central Pacific, (c) east Pacific, (d) south Pacific, (e) 
Indian Ocean, and (f) Atlantic Ocean. Note that the abscissa 
scales are the same except for Fig. 1 0f 

Fig.il. Same as Fig. 10 except for February 1999. 

Fig.12. Contoured-frequency-by-altitude diagrams (CFADs) of monthly 
(February 1998) total Q 2p profiles at 0.5° resolution derived from 
the SLH2 algorithm for various locations. The geographic areas 
are the (a) western Pacific, (b) central Pacific, (c) east Pacific, 
(d) south Pacific, (e) Indian Ocean, and (f) Atlantic Ocean. The 
bin size is 0.5 K. CFAD contour interval: 1% for values < 4% 
and 4% above. Values > 20% are shaded. 

Fig. 13 . Same as Fig. 12 except for February 1999. 




Fig. i. Diagram showing the procedure for refining and validating the 
spectral latent heating (SLH) algorithm. Letters denote convec- 
tive and stratiform classification (C/S), precipitation top height 
(PTH), the precipitation rate at the lowest observable level (P s ), 
and the precipitation rate at the melting level (P m ). The “?” 
means to compare and examine the Q 2 profiles reconstructed 
from the SLH algorithm with the heating profiles from model sim- 
ulations. This consistency check is a necessary precondition for 
the application of the algorithm to actual TRMM PR data. 
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Fig. 2 . Eight-day average profiles of the GCE-simulated heat budget 
[(a) LH, and Q 1R (b) without and (c) with horizontal eddy heat 
flux] and moisture budget [(d) NC, and Q 2 (b) without and (c) 
with horizontal eddy heat flux] for the total (solid), convective 
(dashed), stratiform (dotted) and non-precipitating (thin solid) re- 
gions for the TOGA COARE (19-26 December 1992) case. 
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Fig. 3. Eight-day average profiles of the GCE-simulated moisture bud- 
get [(a) NC, and Q 2 (b) without and (c) with horizontal eddy heat 
flux] for the total (solid), convective (dashed), stratiform (dotted) 
and non-precipitating (thin solid) regions for the TOGA COARE 
(11--17 December 1992) case. 
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Fig. 4. Ensemble-mean, GCE-simulated Q 2p profiles, plotted as func- 
tions of precipitation top height (PTH) from convective (a) and 
stratiform (b) regions, and precipitation rates at the melting level 
from anvil regions (c). Contours indicate values of confidence 
interval for the mean at the 95 % level with Student’s-t test. 
Contour interval is 2.0 K hr 1 for convective and 1.0 K It 1 for 
stratiform and anvil. Thresholds of 0.3 mm h _1 are used for the 
precipitation top detection. 
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Fig. 5. Diagram showing the procedure for deriving latent heating pro- 
files using the spectral latent heating (SLH) algorithm. See the 
text for details. 
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Fig. 6. Eight-day averaged profiles of Q 2p reconstructed by the SLH 
algorithm with the COARE lookup table (thick solid line) and Q 2 
simulated by the GCE model (dotted line) for the (a) COARE 
(19-26 Dec 1992) case, (b) GATE (Sep 1-8 1974) case, (c) 
SCSMEX (Jun 2-9 1998) and (d) KWAJEX (Sep 6-13 1999 ), 
respectively. Leftmost panels are for the convective regions, 
center-left panels for the stratiform regions, center-right panels 
for the non-precipitating regions and the rightmost panels are for 
the total regions. Thin solid lines indicate differences between 
the SLH-reconstructed and the GCE-simulated profiles. 
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Fig. ?. The rms error in the horizontal averaged profiles between the 
SLH algorithm-reconstructed Q 2p and the GCE-simulated Q 2 
for the (a) COARE, (b) GATE, (c) SCSMEX, and (d) KWAJEX 

cases. 
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Fig. 8. Heating from diagnostic calculations (Johnson and Ciesielski 
2002) and the SLH2 algorithm using version 6 of the TRMM PR 
data sets for SCSMEX (15 May - 20 June 1998). 
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Fig. 9. Monthly mean rainfall (mm day -1 ) derived from PR2A25 Ver- 
sion 6 for (a) February 1998 and (b) February 1999. The dry- 
ing/moistening profiles will be compared and examined for the 
various geographic locations identified by the boxes. [From 
Shsge et al. (2007)] 
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Fig. 10 . Monthly (February 1998) mean total, convective and strati- 
form Q 2p profiles derived from the SLH2 algorithm for various 
locations. Total heating profiles derived from the CSH algorithm 
are also shown. The geographic areas are the (a) western Pa- 
cific, (b) central Pacific, (c) east Pacific, (d) south Pacific, (e) 
Indian Ocean, and (f) Atlantic Ocean. Note that the abscissa 
scales are the same except for Fig. 1 0f 
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Fig. ii. Same as Fig. 10 except for February 1999. 







- 6 - 4-20 2 4 6 8 10 12 

(K day-’) 


- 6 - 4-20 2 4 6 8 10 12 

(K day-’) 


- 6 - 4-20 2 4 6 8 10 12 

(K day-') 


Fig. 12 . Contoured-frequency-by-altitude diagrams (CFADs) of 
monthly (February 1998) total Q 2p profiles at 0.5° resolu- 
tion derived from the SLH2 algorithm for various locations. The 
geographic areas are the (a) western Pacific, (b) central Pacific, 
(c) east Pacific, (d) south Pacific, (e) Indian Ocean, and (f) 
Atlantic Ocean. The bin size is 0.5 K. CFAD contour interval: 
1% for values < 4% and 4% above. Values > 20% are shaded. 






Height (km) Height (km) 



Fig. 13. Same as Fig. 12 except for February 1999. 






